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Abstract 

We discuss spontaneous symmetry breaking of (l+l)-dimensional (j)'^ theory 
in hght-front held theory using a Tamm-Dancoff truncation. We show that, 
even though hght-front field theory has a simple vacuum state which is an 
eigenstate of the full Hamiltonian, the field can develop a nonzero vacuum 
expectation value. This occurs because the zero mode of the field must satisfy 
an operator valued constraint equation. In the context of (l+l)-dimensional (j)^ 
theory we present solutions to the constraint equation using a Tamm-Dancoff 
truncation to a finite number of particles and modes. We study the behavior 
of the zero mode as a function of coupling and Fock space truncation. The 
zero mode introduces new interactions into the Hamiltonian which breaks the 
Z2 symmetry of the theory when the coupling is stronger than the critical 
coupling. We investigate the energy spectrum in the symmetric and broken 
phases, show that the theory does not break down in the vicinity of the critical 
coupling, and discuss the connection to perturbation theory. Finally, we study 
the spectrum of the field (p and show that, in the broken phase, the field is 
localized away from = as one would expect from equal-time calculations. 
We explicitly show that tunneling occurs. 



I. Introduction 



Two fundamentally different phenomenological pictures of hadronic matter have de- 
veloped over the past twenty years. One is the constituent quark model where hadrons 
are pictured as being made of a few heavy quarks bound by a confining potential. 
Here, Zweig's rule implies that very little of a hadron's momentum is carried by glu- 
ons. The other is the quark-parton model which is valid at higher energy scales. 
In this picture, a hadron is composed of an indefinite number of light quarks and 
most of the hadron momentum is carried by gluons. We believe that the underlying 
theory for both of these disparate pictures is the non-abelian gauge theory Quantum 
Chromodynamics (QCD). It is believed that the light-front formulation of QCD is 
our best hope of joining the constituent quark model with the parton model and 
QCD. In this elegant approach the field theory is quantized on the null plane surface 
= (x° + x^)/-\/2 = instead of the usual equal-time surface x° = 0. This formula- 
tion avoids many of the difficult problems that appear in the equal-time formulation 
of field theory. 

In the light-front formalism the energy operator does not have a square root op- 
erator, and, consequently, the vacuum structure is relatively simple. Dirac in 
1949, showed that a maximum number of Poincare generators become independent of 
the dynamics in the "front form" formulation, including certain Lorentz boosts. The 
eigenstates of the light-front Hamiltonian have Lorentz scalars = as eigenval- 
ues and describe bound states of arbitrary four-momentum and invariant mass M, 
allowing easy computation of scattering amplitudes and other dynamical quantities. 
However, the most remarkable feature of this formalism is the apparent simplicity 
of the vacuum. The Fock space vacuum is an eigenstate of the full Hamiltonian. In 
other words, the bare vacuum is the physical vacuum. More detail is given in a recent 
review [||. 

For the past several years an increasingly large group of physicists has attempt- 
ing to combine Tamm-Dancoff procedures with light-front quantization in order 
to perform non-perturbative calculations in quantum field theory. The assumptions, 
consistent with the constituent quark model picture, are that a hadron is well ap- 
proximated by only a few particles and that adding more particles only refines this 
initial approximation. This is in stark contrast to the equal-time formulation of QCD 
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where an infinite number of gluons are essential to construct even tlie vacuum 0. If 
these efforts are successful, they could lead to procedures for calculating not only the 
hadron mass spectrum but all the quantities which depend on hadron wavef unctions: 
structure functions, fragmentation functions, et cetera. 

The origin of this remarkably simple vacuum in the light-front formalism is that 
conservation of requires J2i kf = 0. However, > for massive particles, and 
the total light-front Hamiltonian annihilates the perturbative vacuum. In contrast, 
the physical vacuum in equal-time quantization is a highly complex composite of pair 
fluctuations which is believed to produce all phenomena that require long range order: 
color confinement, chiral symmetry breaking, the Goldstone pion, et cetera. How can 
one obtain such non-trivial properties in the light-front formulation of field theory? 
These phenomena occur in light-front field theory because the field includes a zero 
mode operator which is not an independent degree of freedom. This mode is a 
complicated operator-valued function of all other modes in the theory and, since it 
carries zero momentum, it can provide long range order. 

This problem has recently been attacked from several directions. The question 
of whether boundary conditions can be consistently defined in light-front quantiza- 
tion has been discussed by McCartor and Robertson ^ and Lenz 0. They have 
shown that for massive theories the energy and momentum derived from light-front 
quantization are conserved and are equivalent to the energy and momentum one 
would normally write down in an equal-time theory. In the analyses of Lenz et 
al. @ and Hornbostel one traces the fate of the equal time vacuum in the limit 

— s> oo and equivalently in the limit 9 ^ -k /2 when rotating the evolution parame- 
ter T = x'^ cos 9 + sin 9 from the equal-time frame to the light-front frame. Heinzl et 
al. 1^ have considered 0^ theory in (l+l)-dimensions and solved the zero-mode con- 
straint equation by truncating the equation to one particle and retaining all modes. 
They implicitly retain a two particle contribution in order to obtain finite results. 
(We will discuss this in detail later.) They find a critical coupling for the onset of 
spontaneous symmetry breaking. Other authors |T^ find that, for theories allowing 



spontaneous symmetry breaking, there is a degeneracy of light-front vacua and the 
true vacuum state can differ from the perturbative vacuum through the addition of 
zero mode quanta. In addition to these approaches there are many others |jTl], [1^, p^ . 
We present here a discussion of the zero-mode constraint equation for (1+1)- 
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dimensional 0'* field theory with symmetric boundary conditions and show 

how spontaneous symmetry breaking occurs within the context of this model. This 
work builds upon the investigation of the one mode case by Bender, Pinsky, and van 
de Sande in Ref. |T^. Our basic approach is to apply a Tamm-Dancoff truncation 
to the Fock space. We restrict the number of particles and the number of modes to 
be less than some finite limit. This means that we can represent any operator by 
a finite matrix and solve the operator valued constraint equation numerically. The 
truncation assumes that states with a large number of particles or large momentum 
do not have an important contribution to the zero mode. 

We find the following general behavior: for small coupling (large where 
g oc 1/coupling) the constraint equation has a single solution and the field has no 
vacuum expectation value (VEV). As we increase the coupling (decrease g) to what 
we will call the "critical coupling" (^critical, two additional solutions which give the 
field a nonzero VEV appear. These solutions differ only infinitesimally from the first 
solution near the critical coupling, indicating the presence of a second order phase 
transition. Above the critical coupling {g < (^critical), there are three solutions: one 
with zero VEV, the "unbroken phase," and two with nonzero VEV, the "broken 
phase." As shown in Fig. 1, one can plot the VEV as a function of g. We will call 
these the "critical curves." 

In the weak coupling limit {g large) the solution to the constraint equation can be 
obtained in perturbation theory. This solution does not break the Z2 symmetry and 
is believed to simply insert the missing zero momentum contributions into internal 
propogators. This must happen if light-front perturbation theory is to agree with 
equal-time perturbation theory ITsf. This has been shown for 0^ theory in (1+1)- 
dimensions ]16| and we believe that it happens here as well. 

Since the vacuum in this theory is trivial, all of the long range properties must 
occur in the operator structure of the Hamiltonian. Above the critical coupling [g < 
5'criticai) quantum oscillations spontaneously break the Z2 symmetry of the theory. In a 
loose analogy with a symmetric double well potential, we have two new Hamiltonians 
for the broken phase, each producing states localized in one of the wells. The structure 
of the two Hamiltonians is determined from the broken phase solutions of the zero 
mode constraint equation. We find that the two Hamiltonians have equivalent spectra. 
In a discrete theory without zero modes it is well known that, if one increases the 
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coupling sufficiently, quantum correction will generate tachyons causing the theory to 
break down near the critical coupling. Here the zero mode generates new interactions 
that prevent tachyons from developing. In effect what happens is that, while quantum 
corrections attempt to drive the mass negative, they also change the vacuum energy 
through the zero mode and the diving mass eigenvalue can never catch the vacuum 
eigenvalue. Thus, tachyons never appear in our spectra. 

Another way to investigate the zero mode is to study the spectrum of the field 
operator cj). Here we find a picture that agrees with the symmetric double well po- 
tential analogy. In the broken phase, the field is localized in one of the minima of the 
potential and there is tunneling to the other minimum. 

In Sec. II we discuss the derivation of the zero-mode constraint equation and the 
formal operator constraint equation that relates the zero mode to all the other modes 
in the problem. In Sec. HI we apply a truncation allowing only the first nonzero mode 
and many particles. In Sec. IV we examine the truncation where we allow many modes 
and only one particle. Next, in Sec. V, truncations with several particles and several 
modes are examined. In Sec. VI, we examine the spectrum of the field operator 
-\/4vr0. Finally, in Sec. VII we discuss our results and the remaining work needed on 
this problem. 



II. Quantization 

The details of the Dirac-Bergmann prescription and its application to the system 



considered in this paper are discussed elsewhere in the literature |T^, [T^. In this 
section, we will summarize those results and introduce our notation. We define light- 
front coordinates = ± x^) /V2. For a classical field the (0^)^^^ Lagrangian 
is 2 

£ = - - . (2.1) 

We put the system in a box of length d and impose periodic boundary conditions. 
For most of our discussion we work in momentum space. We define qk by 

0(x) = i=E9n(x+)e^'="-- , (2.2) 

V« n 

where = 27Tn/d and summations run over all integers unless otherwise noted. 
Following the Dirac-Bergmann prescription, we can identify first-class constraints 
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which define the conjugate momenta and a secondary constraint which determines 
the "zero-mode" go in terms of the other modes in the theory. This result can also be 
obtained by integrating the equations of motion in position space or differentiating 
the Hamiltonian with respect to the zero mode IP . 



Quantizing, we define creation and annihilation operators a|. and by, 



Atc \k 



afc , afc = aLfc , k 



(2.3) 



which satisfy the usual commutation relations 



[ofc, ai] = 



O-k, 0-/ 



JkA 



kj>0 . 



(2.4) 



Likewise, we define the zero mode operator 




An 



It is useful to define the quantity. 



1 



E 



-'n+«2H hin,0" 



il,i2,...,in^0 



«1«2 ■ ■ - ^r 



(2.5) 



(2.6) 



General arguments suggest that the Hamiltonian should be symmetric ordered 
However, it is not clear how one should treat the zero mode since it is not a dynamical 
field. As an ansatz we will treat Oq as an ordinary field operator when symmetric 
ordering the Hamiltonian. 



Rescaling P , the quantum Hamiltonian is 



Xd 

9^2 , % 



|a^ + ^ + (7S2 + 6S4 



\ — ^ / 2 2 2 



^ k,i,mj^o \j\klm\ 

- c. 



(2.7) 
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where g = 247r/i^/ A. We have removed tadpoles from the symmetric ordered Hamil- 
tonian by normal ordering the third and fourth terms and subtracting, 



3 2 \ ^ 1 
- On > -— . 



(2i 



In addition, we have subtracted a constant C so that the VEV of H is zero. Note 
that this renormalization prescription is equivalent to a conventional mass renormal- 
ization and does not introduce any new operators (aside from the constant) into the 
Hamiltonian. The constraint equation for the zero-mode can be obtained by taking a 
derivative of P~ with respect to Oq. Consequently, we symmetric order the constraint 
equation: 

If 3ao' 
\n\ \ 2 



Q = gaQ + al + ^' 



+ 6S3 



(2.9) 



Using the constraint equation, we can rewrite H as: 
H = ^S2 + 6E4-^ + ^^^(a„a2 
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klm\ 



c . 



(2.10) 



It is clear from the general structure of ( |2.9| ) that (2q clS 3i function of the other modes 
is not necessarily odd under the transform —at, k ^ associated with the 

Z2 symmetry of the system. Consequently, the zero mode can induce Z2 symmetry 
breaking in the Hamiltonian. Along with the zero mode and the Hamiltonian, E„ 
commutes with the longitudinal momentum operator, 











. 



(2.11) 



In order to render the problem tractable, we impose a Tamm-Dancoff truncation 
on the Fock space. Define M to be the number of nonzero modes and N to be the 
maximum number of allowed particles. Thus, each state in the truncated Fock space 
can be represented by a vector of length S = (M + N)\/ {MINI) and operators can 
be represented hj S x S matrices. One can define the usual Fock space basis. 



Fi,n2, 



'fill 



'n2\ 



|0) 



(2.12) 
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where ni + n2 + . . . + um < N . Thus, for example, 

2n 

P^\ni,n2,...,nM) = — (ni + 2n2 H h Mtim) rig, Hm) (2.13) 

E2 |ni,n2, . . . ,nM) = ( y Y ^ ^ ) I'^i, ^^2, • • • , ^m) • (2.14) 

In matrix form, oq is real and symmetric. Moreover, it is block diagonal in states of 
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equal eigenvalue. Due to the Fock space truncation, the commutator 
somewhat different than its full Fock space counterpart. We will use the full Fock 
space commutators to move creation operators to the left as much as possible in a 
product before performing any matrix multiplications. 

III. One Mode, Many Particles 



In this section we review the results of Ref. |14|. Consider the case of one mode M = 1 



and many particles. In this case, the zero-mode is diagonal and can be written as 

N 

ao = fo\0){0\ + Y.h\k){k\ • (3.1) 

k=l 

Note that in ( ^.Ij ) is even under a^. —>■ —a^, k ^ and any non-zero solution breaks 
the Z2 symmetry of the original Hamiltonian. The VEV is given by 

(01010) = ^(0|ao|0) = ^/o . (3.2) 



Substituting ( p.l|) into the constraint equation (|2.9| ) and sandwiching the constraint 
equation between Fock states, we get a recursion relation for {/„}: 

= 9fn + fn + (4n - l)fn + iji + I) + n/„„i (3.3) 

where n < N, and we define /n+i to be unknown. Consequently, {/i, /2, . . . , /n+i} 
is uniquely determined by a given choice of g and /q. In particular, if /q = all the 
/fc's are zero independent of g. This is the unbroken phase. 

Our first objective is to determine which solutions of the recursion relation are 
consistent with the Tamm-Dancoff truncation. If /„ is large for large n then the 
high energy levels will be strongly affected. The paradigm for spontaneous symmetry 
breaking is the symmetric double well potential in ordinary quantum mechanics which 
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has a ground state centered in either well rather than at the symmetry point. This 
paradigm indicates that the behavior of the system is unaffected by the barrier for 
energies far above the barrier separating the wells. Hence, we only seek solutions 
where /„ is small for large n. We caution here that this paradigm does not translate 
entirely to the light-front formulation since the symmetry breaking occurs in the 
Hamiltonian and the vacuum itself is unaffected. 

We begin our analysis of Eq. ( |3.3| ) by dropping the cubic term and finding solutions 
to the resulting linear system. Assume for now that N = oo and that is decreasing 
with n. For large n, the terms linear in n dominate, Eq. (|3.3|) becomes 

/n+i+4/„ + /„_i = 0. (3.4) 

There are two solutions to this equation: 

/„oc(v^±2)". (3.5) 



We reject the plus solution because it grows with n. Dropping the cubic term from 
(|37^ ) we define the generating function 



Hz) = E /n^" • (3-6) 

n=0 

If fn goes like (^y/S — 2^ then the radius of convergence of F{z) is 2 + -\/3 and we 
expect F{z) to be singular at |2;| = 2 + Similarly, if /„ ~ ("\/3 + 2^ , then we 
expect F{z) to be singular at |2;| = 2 — 

The function F{z) satisfies a differential equation whose solution is 



F{z) fz + 2-^/3\ fz + 2 + V^\ 



(3.7) 



F{0) \ J \ 2 + J 

Note that this solution for F{z) has singularities at the expected values of z. If 
we want fn to have the asymptotic behavior (^^/S — 2^ for large n, then we must 
eliminate the branch point of F(z) at \z\ = 2 — v^- This gives the condition 

- '^~^~'' =K, K = 0,l,2... (3.8) 
Concentrating on the K = case, we find a critical coupling 

^'critical = 3 - \/3 (3.9) 
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or 

Acriticai = 47r (3 + v^) /i^ ~ 60/i2, (3.10) 

consistent with equal-time calculations p9|] . 

The solution to the linearized equation is an approximate solution to the full 
equation ( |373| ) for /o sufficiently small. We need to determine solutions of the full 
nonlinear equation which converge for large n. We will use both the 5-expansion and 
numerical methods to do this. 



The (5-expansion is a powerful perturbative technique for linearizing nonlinear 
problems. It has been shown to be an accurate technique for solving problems in 
differential equations, quantum mechanics, and quantum field theory pO| . 

We rewrite Eq. ( |3.3|) as 

{g-l+ An) U + fn^'' + (n + 1) U+i + = . (3.11) 

Setting 6 = gives the linear finite difference equation which is the zeroth-order 
approximation in the 5-expansion. One then expands in powers of 6 about 6 = 0. 
One recovers the problem of interest aX 6 = 1. 

We find to first order the critical curve for 5 = 1, 

g=(2-V3) (^l + -^ln(2 + v^)^ -ln/2. (3.12) 

These results are plotted as the dashed curve in Fig. 1. The expansion behaves badly 
near fo = because of the Iu/q term in the expansion. The ^-expansion analysis 
clearly shows that there is a critical curve and not merely a critical point. 



We can also study the critical curves by looking for numerical solutions to Eq. (|3.3|). 
The method used here is to find values of /o and g such that /tv+i = 0. Since we seek 
a solution where /„ is decreasing with n, this is a good approximation. We find that 
for g > 3 — y/S the only real solution is /„ = for all n. For g less than 3 — \/3 there 
are two additional solutions. Near the critical point |/o| is small and 

/„^/o(2-v^)" . (3.13) 

The critical curves are indicated by the solid lines in Fig. 1. These solutions con- 
verge quite rapidly with N. The critical curve for the broken phase is approximately 
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parabolic in shape: 

^ ^ 3 - V3 - 0.9177/o^ . (3.14) 



It is instructive to study the behavior of the constraint equation (|3.3| ) away from 
the critical curves. In Fig. 2 we plot as a function of n and /o for g = 1.2. We 
see that, as n becomes large, all the increase and as /o approaches the critical 
curve, which is at /o ~ 0.2700 for g = 1.2, all the |/n|'s decrease rapidly. As /o 
increases beyond the critical curve the l/nl's increase rapidly once again. The fact 
that \fn\ increases rapidly on both sides of the critical curve is a manifestation of the 
nonlinearity in ( p.3| ). 

In Fig. 3 we extend the critical curves to the negative g region. This region corre- 
sponds to negative values of /i^ in the Hamiltonian. We have not studied the spectrum 
for fi"^ < {g < 0); however, it is easy to see from the figures that nothing unusual 
happens, at least for small g. As predicted by Eq. ( |3.8| ), there are additional solutions 
near the critical curve for the unbroken phase. The curves shown are independent of 

for A^ large. It is not clear how to interpret these solutions. 

We can also study the eigenvalues of the Hamiltonian for the one mode case. The 
Hamiltonian is diagonal for this Fock space truncation and, 

(-1 H \n) = ln{n - 1) + - | - + + ^fl, - C . (3.15) 

The invarient mass eigenvalues are given by 

I / I / 247r ' ^ ^ ' 

In Fig. 4 the dashed lines show the first few eigenvalues as a function of g without the 
zero-mode. When we include the broken phase of the zero mode, the energy levels 
shift as shown by the solid curves. For g < ^fcriticai the energy levels increase above 
the value they had without the zero mode. The higher levels change very little, as 
our paradigm would suggest, because /„ is small for large n. 

IV. Many Modes, One Particle 

Now, let us look at the case of one particle A^ = 1 and many modes. In this case, the 
zero mode is diagonal and can be written as 

M 

ao = bo\0) {0\ + Y.hai\0) {0\a, . (4.1) 

k=l 
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Substituting this into the constraint equation ( |2.9| ) and sandwiching the constraint 
equation between various states, one obtains a system of equations for the coefficients 

= hl + gh, + Y.-^ (4.2) 

fc=i 

The last term of ( |4.3| ) couples the one and two particle sectors of a^. Assuming that 
the Fock space truncation is a good approximation, we set this term equal to zero. 
Note that 6^ = 0, = 0, 1, . . . , M, is always a solution of this system of equations; 
this is the unbroken phase. Evaluating the system of equations numerically, we find 
the critical curves shown in Fig. 5. 

Let us calculate the critical coupling. As with the one mode case, we drop the 
cubic terms from ([4.2|) and ([4.3|) and find a solution to the resulting linear system. 
Defining 9 = g — Y^kLi ^/k, we find 

bo 



md + 4 



(4.4) 



M 



e=y] —— — - (4.5) 

or, in the limit of large M, 

4^ = 7^ + ^(l + 4/^)+0(l/M) = 0.6267537... + 0(1/M) (4.6) 

where 7^ is Euler's constant, and \l/(x) is the digamma function. Thus, the critical 
coupling is logarithmically divergent, 

^critical = logM + 7^ + + 0(1/M) . (4.7) 

We compare this expression with numerical results in Fig. 6. 

Next we need to examine whether the Fock space truncation is consistent, i. e. 
whether bm converges for large m. If we look at the linearized solutions, Eq. (^4.4|), 



we see that does converge near the critical coupling. Now, examine the full non- 
linear equations. From Eq. (|4.3| ), we see that bm is independent of m if 6^ = —bo/ 4:. 
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Substituting this into Eqs. ( |4.2| ) and ( [4 .31 ), one can solve for g. Define to be this 
value of g, 

k=l 

Note that g^ < (^critical for all M and that go is logarithmically divergent in M. For 
g = go, bm is independent of m. For g < go, bm diverges with m, and, for g > go, 
bm converges with m. This means that the solution to the constraint equation is 
inconsistent with the Fock space truncation for g < go- 

Elsewhere it has been suggested that one use a somewhat different ansatz for the 
zero mode 0, 

M 

ao = Co + Yl Cfealflfc • (4.9) 

k=l 

How does it compare with the definition in Eq. ( |4.1| )? Let us look at some matrix 
elements in each case: 

(0|ao|0) = 6o (0|ao|0) = co 

(0| akaoal |0) = bk (0| akdoal |0) = cq + Ck 

(0| a^aiaoalal |0) = (0| akaidoolal |0) = Cq + + Q , k ^ I 

^ (0| (a.fao (4)' |0) = i (0| {a^f do (4)' |0) = Cq + 2ck . (4.10) 

The important point is that this ansatz makes very different assumptions about the 
two particle sector. If we identify bo with cq and bk with co + Ck, it would be equivalent 
to a different choice choice for the last term in equation ( [4.3| ): 

(0| amakaoaWm |0) _ y> ^fc - ^ ^ y~ 1 ^ '^bm - bo 

fe=l *^ k=l k=l'^ ^ 

Using this assumption for the two particle sector gives us a different system of equa- 
tions for the coefficients, 

= bl + gbo + Y.-^ (4.12) 

fc=i 

= bt + gh^ + -h^ + Y.-^- (4.13) 

These are equivalent to equations (3.6a) and (3.6b) in Heinzl, et al. [§]. In this case, 
the zero mode matrix elements along with the critical coupling converge in the limit 
of large M. 
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How do we interpret this ansatz in terms of a Fock space truncation? We have 
assumed, in essence, that the zero mode has a contribution from the 2 particle sector 
that cannot be neglected. Is this the correct way to truncate the Fock space? We 
will answer this question in Sec. V. 

We turn our attention to the Hamiltonian. For this Fock space truncation the 
Hamiltonian is diagonal. Substituting Eq. ( [4.1| ) into Eq. ( |2.10| ) and sandwiching 
between states, we obtain 

(OW) = {H-^-i-C (4.14) 

k=l 

1,4 1,2 Ml 7 2 

{^\a^Ham = ^-^--^-blY^Z + J^-C- (4-15) 
m 4 2m ^ k Am 

We choose C such that the VEV of H is zero and plot a spectrum in Fig. 7. Note 
the vertical line representing the value of go- Presumably, to the left of this line, the 
Fock space truncation breaks down. Also, note that states with larger energy have 
smaller longitudinal momentum P+. Clearly, the zero mode has a large effect on the 
spectrum. 



V. Many Modes and Many Particles 

We turn our attention to the more general case of many modes and many particles. 
Many of the features that were seen in the one mode and one particle cases remain. 
Plotting the VEV vs. g, one finds that the critical curves have the same form with an 
unbroken phase, a broken phase, and a critical coupling. There is rapid convergence 
as one increases and a logarithmic divergence as one increases M. However, now 
one can have more than one state of a given momentum. Thus, H and ao are no 
longer strictly diagonal and S3 is nonzero. 

Consequently, the zero mode for the unbroken phase is no longer zero. We can use 
Eq. ( p.9|) and solve for the unbroken phase of perturbatively in A or, equivalently, 
in 1/g: 

ao = --E3 + f 2S,S3 + 2S3S, + E »^^3aI + 4S3a.-S3 \ ^ ^^^^^3^ _ 
9 9 \ k=i J 
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For each order in one can see that is odd under the transform at — > — ctfc, 
A; 7^ 0. Consequently, 

(a|ao|/3) = if (a|iV|a) - {I3\N\I3) is even (5.2) 

where is the number operator and \a) and \(3) are Fock space basis states. In 
particular, the diagonal matrix elements of ao are zero. It is important to note that, as 



expected, this expansion does not produce the broken phase ||10|]. When substituted 
into the Hamiltonian, the perturbative expansion of the zero mode produces new 
interactions in perturbation theory. It is generally believed that equal-time and light- 
front perturbation theories are equivalent |jl5| and that the zero modes will change 



the internal propagator to be the full propagator that one would have if there were 
a dynamical zero mode. Maeno has, in fact, shown that this is true for (jr" theory in 



l+l)-dimensions 16 



In order to calculate the zero mode for a given value of g one converts the constraint 
equation ( p.9| ) into an 5 x 5 matrix equation in the truncated Fock space. This 
becomes a set of S"^ coupled cubic equations and one can solve for the matrix elements 
of ao numerically |21[]. Considerable simplification occurs because ao is symmetric 
and is block diagonal in states of equal momentum. For example, in the case M = 3, 

= 3, the number of coupled equations is 34 instead of S"^ = 400. In order to 
find the critical coupling, we take (0|ao|0) as given and g as unknown and solve the 
constraint equation for g and the other matrix elements of ao in the limit of small but 
nonzero (0|ao|0). Critical coupling as a function of M and is shown in Fig. 8. We 
see that the solution seems to be consistent with our earlier results: there is quick 
convergence as N increases and a logarithmic divergence as M increases. 

Thus, we believe that the ansatz introduced in Eq. ( [4.9|) produces misleading 
results. If one were to extend the ansatz to include more particles, the convergent 
result for the critical coupling would disappear. 

As an example, let us examine some matrix elements of ao for the truncation 
M = 3 and A^ = 3. In this case, 5* = 20 and we restrict our attention to states with 
P+ < 3(27r)/(i (numerical calculations are performed with all 20 states). Arrange the 
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basis states in the following order: 





=0 


|0,0,0) 


P+-- 


~^ d 


|1,0,0) 


P+-- 




|0,1,0) 






|2,0,0) 


P+-- 




|0,0,1) 












|3,0,0) 



(5.3) 



Given g — 2.8 above the critical coupling gicriticai = 2.301, we solve the constraint 
equation, 

/ 







-0.287 
-0.287 



-0.281 
-0.281 -0.458 
-0.458 



V 



(5.4) 



which leads to a Z2 symmetric Hamiltonian. Now, we lower the coupling to = 1.8. 
The solution for the unbroken phase is now 



/ 



-0.318 
-0.318 



-0.315 

-0.315 -0.499 
-0.499 



V 



(5.5) 



In addition, we have two solutions for the broken phase. There is a solution with 
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(0|ao|0) > 0, 

/ 0.787 



ao 



-0.246 



-0.251 
-0.314 



-0.314 
0.072 



-0.270 
-0.310 
0.005 



-0.310 
0.070 
-0.499 



0.005 
-0.499 
-0.019 



V 

along with the (0|ao|0) < solution, 
/ -0.787 



(5.6) 



ao 



0.246 



0.251 -0.314 
-0.314 -0.072 



0.270 

-0.310 
-0.005 



-0.310 
-0.070 
-0.499 



-0.005 
-0.499 
0.019 



V 



(5.7) 



When we substitute the solutions for the broken phase of ao into the Hamiltonian 
( 2.10| ) we get two Hamiltonians and corresponding to the two signs of (0|ao|0) 
and the two branches of the curve in Fig. 1. This is the new paradigm for spontaneous 
symmetry breaking: multiple vacua are replaced by multiple Hamiltonians. Picking 
the Hamiltonian defines the theory in the same sense that picking the vacuum defines 
the theory in the equal-time paradigm. The two solutions for ao are related to each 
other in a very specific way. Let H be the unitary operator associated with the 
Z2 symmetry of the system; Ha^H^ = — a^, A; 7^ 0. We break up ao into an even 
part Ha^H^ = and an odd part Ha^H^ = — a^. The even part af breaks the 
Z2 symmetry of the theory. For g < (/critical, the three solutions of the constraint 
equation are: corresponding to the unbroken phase, a^ + af corresponding to 
the (0|ao|0) > solution, and — for the (0|ao|0) < solution. Thus, the two 
Hamiltonians are 

H+ = H[ak,a'^ + a^) (5.8) 

and 

H^a^.a'^-al) (5.9) 



H 
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where H has the property 



H {ak, ao) = H {-ak, -ao) 



(5.10) 



and ak represents the nonzero modes. Since 11 is a unitary operator, if is an eigen- 
vector of H with eigenvalue E then n|\l/) is an eigenvalue of UHlV with eigenvalue 



E. Since, 



= H{ak,a'^ + a^) = , 



O E 

-ak, — Oq — Oq 



(5.11) 



and H~ have the same eigenvalues. 

Using the M = 3, = 3 case as an example, let us examine the spectrum of H. 
Removing the zero mode entirely, the rescaled Hamiltonian, Eq. ( |2.10| ) is 



/ 



H 



\ 



9 
2 







3 + 2g 



3. 

3 





3 + 
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V2 



V2 


9 + 3g 



V 



(5.12) 



For large g the eigenvalues are obviously: 0, g, g/2, 2g, g/3, 3g/2 and 3g. However as 
we decrease g one of the last three eigenvalues will be driven negative. This signals 
the breakdown of the theory near the critical coupling when the zero mode is not 
included. 

Including the zero mode fixes this problem. Fig. 9 shows the spectrum for the three 
lowest nonzero momentum sectors. This spectrum illustrates several characteristics 
which seem to hold generally (at least for truncations we have examined, N + M < 6). 
For the broken phase, the vacuum is the lowest energy state, there are no level 
crossings as a function of g, and the theory does not break down in the vicinity 
of the critical point. None of these are true for the spectrum with the zero mode 
removed or for the unbroken phase below the critical coupling. The lowest eigenvalue 
has a minimum at or near the critical point and the minimum appears to decrease 
with the number of modes and particles. A more precise statement will have to await 
a fully renormalized treatment. 
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VI. Spectrum of the Field Operator 



So far wc have examined the zero mode and its effect on the spectrum of the Hamilto- 
nian. How does the zero mode affect the field itself? Since is a Hermitian operator 
it is an observable of the system and one can measure for a given state \a). The 
result is given by simple quantum mechanics. Let us define the eigenvalues 0j and 
eigenvectors \xi) of y/Ancf) : 

A/47r0 \Xi) = 4>i\Xi) , {Xi\Xj) = ■ (6.1) 

The probability of obtaining 0j as the result of a measurement of y/^ncp for the state 

I") is \{Xi\a)f- 

In the limit of large N. the probability distribution becomes continuous. If we 
ignore the zero mode, the probability of obtaining as the result of a measurement 
of -\/47r0 for the vacuum state is 

where r = J2kLi ^/k. The probability distribution comes from the ground state wave 
function of the Harmonic oscillator where we identify with the position operator. 
This is just the Gaussian fluctuation of a free field. Note that the width of the 
gaussian diverges logarithmically in M. When is finite, the distribution becomes 
discrete as shown in Fig. 10. 

In general, there are + 1 eigenvalues such that (Xi|0) 7^ 0, independent of M . 
Thus if we want to examine the spectrum of the field operator for the vacuum state, 
it is better to choose Fock space truncations where N is large. With this in mind, we 
examine the = 50 and M = 1 case as a function of g in Fig. 11. Note that near the 
critical point. Fig. 11a, the distribution is approximately equal to the free field case 
shown in Fig. 10. There is no symmetry breaking and the field is symmetric about 
zero. As we move away from the critical point. Figs, llb-d, the distribution becomes 
increasingly narrow with a peak located at the VEV of what would be the minimum 
of the symmetric double well potential in the equal-time paradigm. In addition, there 
is a small peak corresponding to minus the VEV. In the language of the equal-time 
paradigm, there is tunneling between the two minima of the potential. 
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VII. Discussion 



In the context of (0"^) 



1+1 



on the light front, the vacuum of the full theory is always the 



perturbative vacuum. The zero mode, which satisfies an operator valued constraint 
equation, produces the long range physics of the theory including spontaneous break- 
ing of the Z2 symmetry. We have found that the constraint equation can be solved 
using a Tamm-Dancoff truncation of the Fock space. Even for the one mode trun- 
cation, we find a critical coupling consistent with the best equal time calculations. 
Increasing the Tamm-Dancoff truncation, we find rapid convergence with the total 
number of particles N . This is to be contrasted with the equal-time approach where 
an infinite number of particles are required to produce a critical point. In the weak 
coupling limit, the zero mode ensures that light-front perturbation theory agrees with 
equal-time perturbation theory. Above the critical coupling the zero mode develops a 
contribution that breaks the symmetry of the theory. There are two such solutions 
to the constraint equation: one with (O|0|O) > and one with (O|0|O) < 0. This, in 
turn, gives rise to two Hamiltonians which have the same spectrum. The spectrum 
has the expected behavior: the Fock vacuum is the state of lowest energy and the 
lowest eigenvalue has a minimum at the critical coupling. Closer inspection of the 
field shows that tunneling occurs between positive and negative eigenvalues in the 
broken phase. 

In this work, we apply a simple mass renormalization to the symmetric ordered 
Hamiltonian. This prescription removes tadpoles from ordinary interaction terms 
and would properly renormalize the theory if the zero mode were removed. However, 
inclusion of the zero mode produces logarithmic divergences in the constraint equation 
and in the resulting Hamiltonian. This can be seen from Figs. 6 and 8 where (^critical 
grows with the number of modes. The number of modes is equivalent to a large 

cutoff in a standard discretized light-front quantization calculation. We therefore 
believe that this behavior is a renormalization effect related to the fact that the new 
interactions are not normal ordered. To a very good accuracy, 



This type of logarithmic growth indicates that the system needs further renormaliza- 
tion. The most direct approach would be to find some non-perturbative method of 




(7.1) 
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normal ordering the resulting Hamiltonian. Another approach is to add a constant to 
g and multiply A by a constant so that logarithmic divergences are removed. In the 
one particle case, this will work. However, it is unclear if we can successfully extend 
it to the more general case > 1. Another possibility is to add more operators to 
the Hamiltonian. For instance, if we add a term of the form / dx~ dy~(f){x)(f){y), we 
can remove the logarithmic divergence in the constraint equation. 

For {(f)'^)l_^_l we believe that a vanishing of the mass gap is associated with the 
critical couphng 0. In Fig. 7 we see that the gap between the vacuum and the 
lowest energy excited state is minimized at the critical coupling. One could imagine 
that, in the limit of large M, the gap between the vacuum and the first excited state 
goes to zero at the critical coupling. However, the first excited state is the one with 
the largest longitudinal momentum and is always at the "edge" of our Fock space 
truncation in M. Thus, our Fock space truncation may not be particularly well 
suited for an investigation of the vanishing of the mass gap. 

We chose as a model for its simplicity. One problem with this choice 

is that the scalar field is dimensionless. Power counting arguments suggest that 
arbitrary powers of are allowed in the Hamiltonian. Consequently, the zero mode 
can be a very complicated object for this model. In contrast, for 0^ theory in (3+1)- 
dimensions the only allowed local operators are 0", n < A. For physically interesting 
theories like QCD, the operators that are allowed by power counting arguments are 



much more restricted However, they still allow a very large set of operators and 
are therefore in some sense similar to the (0^)^^^ example. Solution of the zero mode 
problem for the model we chose may be conceptually more difficult. 

In theories like QCD we expect to have, in addition, zero modes that are dynamical 
degrees of freedom. A recent study of pure glue QCD in (l+l)-dimensions shows that 



the zero mode of A'^ has this property |T3|, Thus, in QCD, we expect to have 



dynamical as well as nondynamical zero modes making the problem quite complicated. 
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Figure Captions 



Figure 1. /o = v47r(O|0|O) vs. g = 247r/i^/A in the one mode case. The sohd 
curves are the critical curves obtained from numerical solution of Eq. ( |3.3| ) with 
= 10. The dashed curve is the critical curve obtained from the first-order 
(5-expansion. 

Figure 2. vs. n and /o for g = 1.2 using Eq. ( p.3|) . 

Figure 3. Critical curves for a larger range of g in the one mode case with = 10. 
The curves shown are independent of N, N > 10. 

Figure 4. The lowest three energy eigenvalues for the one mode function 
of g from the numerical solution of Eq. ( p.l5|) with = 10. The dashed lines 



are for the unbroken phase /o = and the solid lines are for the broken phase 



Figure 5. bo = v47r(O|0|O) vs. g for the one particle case. The critical curves are 
obtained from numerical solution of Eqs. ( |4.2| ) and ( [4.3|) with M = 10. The 
vertical line represents g^. Note the presence of additional solutions for g small. 



Figure 6. (^critical vs. M in the one particle case. The solid curve is from Eq. 

and the points are from numerical solution of (|4.2|) and ( [4.3|) . The dashed curve 
is go from 



Figure 7. Eigenvalues as a function of g in the one particle case using Eq. ( |4.15| ) 



with M = 10. The solid lines are for the broken phase and the dashed lines are 
for the unbroken phase. The vertical line represents ^fo- 

Figure 8. Critical coupling vs. M and N. 

Figure 9. The spectrum for M = 3, = 3 and (a) P+ = l(27r) /d, (b) P+ = 2(27r) /d, 
and (c) = 3(27r)/(i. The dashed line shows the spectrum with no zero mode. 
The dotted line is the unbroken phase and the solid line is the broken phase. 



Figure 10. Probability distribution of eigenvalues of v47r0 for the vacuum with 
M = 1, A^ = 10, and no zero mode. Also shown is the infinite A^ limit from 
Eq. (pD. 
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Figure 11. Probability distribution of eigenvalues of -\/47r0 for the vacuum with 
couplings (a) g ^ 1, (h) g ^ 0, (c) g = -1, and (d) g ^ -2. M ^ 1, N ^ 50, 
and the positive VEV solution to the constraint equation is used. 
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